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Abstract 

Problems involving disordered systems are usually analyzed for systems with 
random disorder. However, there are many systems in which the main disorder 
involves clusters with correlated differences between their properties and those 
of the average system. A new approximation, the average trace approximation, 
is proposed for calculating the diagonal elements of the Green function, and 
hence the density of states, in such systems. As an example, application of the 
method to a simple cubic array of harmonic oscillators shows that correlation 
in the disorder leads to a peak in the low frequency density of states, a result 
confirmed by computer simulations. 
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There are many problems in the physics of condensed matter that involve 
the behavior of disordered systems. The simplest type of disorder involves fluc- 
tuations in the diagonal elements of the system's Hamiltonian. For electronic 
systems, such problems correspond to fluctuations in the electron's self energy, 
while for vibrational problems it corresponds to isotopic impurities, i.e. fluc- 
tuations in the masses of the atoms. The most powerful method to study the 
effect of such disorder is the coherent potential approximation (CPA). 1 Another 
type of disordered system that has been widely studied is that of off-diagonal 
disorder in the elements of the Hamiltonian between adjacent sites. Exam- 
ples of this for electrons include random hopping rates between different sited^ 
and the migration of localized electronic excitations among guest molecules in 
a hostPI The corresponding problem for vibrations is that of bond disorder, 
which is modelled by systems of harmonic oscillators with random distributions 
of spring constants^ All these authors use a coherent medium approximation 
(CMA) that is essentially equivalent to the CPA. Spin glasses are an example 
of a more complicated random system, which we will not consider in this paper, 
since they are not usually treated in terms of Green functions. 

In all the above system, it is usually assumed that the disorder is random. 
Even when correlated disorder is considered, it is often treated in terms of a 
probability distribution of the correlation decaying with distance between the 
sites.^ This is not the type of correlation that is relevant to specific defect clus- 
ters involving several sites or bonds with correlated properties different from 
those of the surrounding medium. Such defects can easily occur, for instance 
with the bonds to atoms with irregular coordination numbers and the electronic 
states associated with them. Other examples include clusters of molecules in 
liquids, and atoms or electronic states in the region of grain boundaries in crys- 
tals or in regions of increased or decreased density in glasses and fluids. For 
electronic systems, the homomorphic coherent cluster approximation (HCPA^l 
was proposed (which actually was the original basis for the CMA) to treat such 
problems, and had some success. However, there are grave practical difficulties 
in using it in more than one dimension for clusters containing more than two 
sites, since it is difficult to find clusters of identical shape (homomorphic) that 
span the whole system and for which calculations are practicable. For instance, 
the smallest such cluster for a simple cubic lattice contains seven sites, the cen- 
tral one and its six nearest neighbors. In addition, a fundamental objection has 
been raised to the use of the HCPA, 7 namely that the sites on the edge of the 
cluster are regarded partly as belonging to the cluster and partly to the average 
surrounding medium, so that it is far from clear when use of the HCPA is justi- 
fied. In this paper, we present a new method for calculating the density of states 
in systems containing a low density of clusters with correlated disorder, which 
is based on the T-matrix formalism for the calculation of Green functions and 
a new approximation to the diagonal elements of the Green function which we 
call the average trace approximation (AVTA). As an example of the application 
of this method, we consider its application to a simple cubic array of harmonic 
oscillators containing defects for which all the spring constants are weak. For 
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this system, we find a low frequency peak in the density of states which is not 
present for the same amount of random disorder, and our results (and hence the 
accuracy and usefulness of the AVTA) are confirmed by computer simulations. 
The systems that we consider have a Hamiltonian of the form 

H = H + H p (1) 

where Hq is a Hamiltonian having a known Green function Go and the pertur- 
bation (not necessarily small) H p can be expressed in the appropriate basis set 
as the sum of Hamiltonians H a , a = 1...M associated with M disjoint clusters 
For convenience, we will consider the most common case, in which Hq is the 
Hamiltonian of a crystal and we use the site representation, in which we denote 
a site by n, so that an operator F can be written as J2 m J2 n |m > F < n|. 
Our aim is to calculate the density of states g(E) in the system from the trace 
of its Green functionP 

g(E) = T-lm{Tr[G(E±ie)}}. (2) 
For the calculation of G, we use the T-matrix formalism, and write 

G(z) = Go(z) + Go(z)T(z)Go(z), (3) 

where 

T{z)=H p [I-G {z)H p }-\ (4) 

For random disorder, the above formalism is usually used either in the average 
T-matrix approximation (ATA), where one writes < G(z) >— Gq(z) + Gq(z) < 
T(z) > Gq(z), or in the CPA, in which one finds the value of z for which 
< T(z) >— 0. However, even the extension of the ATA to clusters of more than 
two sites involves all the above-mentioned problems of the HCPA. In addition, 
even for two sites there are problems in defining the operator P required to 
solve equation (7), as described below after equation (12). In this paper, we 
propose an entirely different method for calculating the diagonal elements of G, 
and hence the density of states. 

While a general operator F has non-zero matrix elements between any pair 
of sites, the Hamiltonian H p of the perturbation is the sum of Hamiltonians 
H a each of which is associated with a cluster containing only a few sites. In 
that case, if a cluster centred around n a contains only the site n Q and some 
adjacent sites n Q + a^, k — 1..R, it is convenient to introduce the row vector 
\a >= (|n Q >, |n Q +ai >, \n a +aR >), as proposed by EconomoiPln order to 
present our method, we consider first the conventional case where each cluster 
consists of only a single site, so that two clusters a = 1,2 contain only |ni > 
and |ri2 > respectively, and write G(2; 111,112) = < ni|G(z)|n2 >, and similarly 
T(z; ni, 112) = < ni|T(z)|n2 >. In this notation, we can write 

M M 

H x = Y,H a = ^ |n Q > V a < n a \. (5) 

a— 1 a — 1 
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where V a is the matrix (in this case a scalar) describing the perturbation at 
site n a . For a system containing only a single defect at site n a , it readily 
follows from equations (4) and (5) that the corresponding operator T a is just 
T a = \n a > t a < n a \, where t a is a scalar, 

t a = V a [l- <n a \Go\n a >V a }-\ (6) 

If there are a number of defects, one can write 

T = J2Q a = UT a + T a G £ Qp], (7) 

a a P=jia 

where Q a is the contribution to the scattering matrix T from the defect at n Q 
in the presence of all the other defects. Iteration of this equation leads to the 
standard equation 

Qa — E "Ta + T a G Tp + T a G a TpGnT~ ( + ... (8) 

These equations can be generalized to defects containing more than one site 
by replacing the scalars |n Q > with the vectors \a > throughout, in which case 
V a , G a =< a\Gn\a > and t a will be square matrices. In order to show this, we 
consider a specific cluster a, and write m*, = n a + a*,. Then, for instance 

H a = \a > V a < a\ = l m fc > V a < m,|, (9) 

k,i 

and since < r|s >= S r ^ s it follows that < r|_ff Q |s >= 5 rymk S Sylni Vki, and so 
that Vu =< mk\H a \mi > is the element (k,l) of the matrix V a =< a\H a \a >. 
Similar definitions apply to the matrices G a =< a\Gn\ot > and t a , and it can 
easily be shown that the equivalent of equation (6) for our matrices is 

t a = V a [I-G a Y a \- 1 (10) 

Thus, in general we can write 

T = V E E "T a + T a G Tp + T a G TpG T, ( + (11) 

where now 

r a = |a>t a <a|=5^t fc j|m fc ><mj|. (12) 

k,l 

However, while for single site clusters equation (7) can be solved by defining 
an operator P = Gq— < n a |Go|n a >, this method does not apply to solving 
equation (11) for clusters containing more than one site. Instead, we note that 
for the density of states we only need the trace of G and so of Go and of GnTGo- 
In order to calculate the latter, we introduce the average trace approximation 
(AVTA), in which we write 

M 

(Tr{G TG }) w (T>{G (]T T a )G }) = M(Tr{G \a > t a < a\G }), (13) 

a=l 
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where (.) denotes the mean value and M is the number of defect clusters. In 
the AVTA it is essential to take the trace before averaging, since (in contrast 
to the ATA and CPA) we do not average out the defects over the whole system 
and calculate (G) for the periodic average system. The justification of the 
AVTA for small defect densities is that the first term which we are ignoring in 
Tr{G TG } is Tr{G T a G TpGo} where ^ a. On taking the trace in the site 
representation, we write 

Tr{G a T a G T p G } = ^ < n|G |a > t a < a\G o \0 > t < 0\G o \n > . 

n 

Each term of this sum involves the matrix elements of Go between sites belonging 
to two different defects multiplied by the matrix elements of Go between a site 
n and one site from each defect. These terms will be small if the defects are far 
apart, as they will be for low defect densities. In terms of Feynman diagrams, 
the first terms that we are ignoring correspond to propagation from an arbitrary 
site to one defect, from that to the second defect, and then back to the original 
site. If we were interested in non-diagonal elements of G, for some of which a 
path starts at a site close to cluster a and ends at a site close to cluster 0, this 
approximation would not be justified. 

For the calculation of Tr{Go|T a |Go} for a given cluster a, we write Gq(z) 
in the form Gq(z) — ^ k |k >< k|/(z — E^) and use equation (12) for T a , while 
< k|n >= exp(— ik • n). It then follows that 

rr sn t„\\T ( \\n (~\\ M S^ST+ exp[ik • (a m - a ; )] 
Tr{Go{z)\T a (z)\G a (z)} = j^}^l^ tlm ( z _ Ek )2 ' <- 14 ) 

l,m k 

which depends only on the distances a m — a; between pairs of sites in a cluster. 
In order to evaluate the trace in terms of the known expressions for Gq(z; aj, a m ), 
we write equation (14) in the form 

Tr{G {z)\T a (z)\G (z)} = -M%- V t lm G (z; a,, a m ). (15) 

oz f—* 

In order to demonstrate the usefulness and accuracy of the AVTA, we now 
apply it to a specific system, and compare its results with those of computer 
simulations. The model system that we consider is a simple cubic lattice of 
coupled harmonic oscillators with spring constants K. For the sake of simplicity 
we chose a binary distribution, with a small fraction of oscillators having a 
spring constant K\ (which we refer to as weak springs) much less than the 
spring constant Kg of most of the oscillators. Extensive results for a random 
distribution of weak springs in this system have been presented by Schirmacher 
and Diezemann. P who found that a concentration of 10% of weak springs had a 
negligible effect on the density of vibrational states. In contrast to their model, 
we considered systems in which the weak springs occur in clusters, with the same 
type of cluster in each defect. The defects that we consider are those in which 
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a site is linked to its six neighbors by weak springs, with Kx/Kq =0.1. The 
results for more general systems of this type, in which a site is linked to its six 
neighbors by z weak springs with spring constant K\ and 6 — z ordinary springs 
with spring constant Ko, as well as the results for other values of Ki/Kq, will 
be reported elsewhere The concentration c = M/N of these defects was chosen 
to be 1.5%, so that the total concentration of weak springs was 9%. Since the 
clusters are identical, it follows from equation (15) that they each give the same 
contribution to Tr{GoT a Go}- Thus it is sufficient to calculate Tr{GoT a Go} for 
one such cluster, so that 

Tr{G} = Tr{G a } + cNTr{G Q T a G a }, (16) 

and for the density of states per site 

^Tr{G} = G (z; 0, 0) - c^- ^ t lm G {z; a /; a m ). (17) 

In order to check the accuracy of the AVTA, we also performed a number of 
simulations on cubes with 13-16 oscillators on each side, and took a smoothed 
average of the density of states obtained for these systems. 

The results of our calculations are shown in figure 1, in which we compare 
the densities of states in the periodic system (no defects) with that in the system 
with a concentration of 1.5% of these defects. We note first that in the periodic 
system the density of states g(uj) is fairly smooth, with the well-known van-Hove 
singularities, where there is a discontinuity in the gradient of g(ui) : at uj = 2 
and uj = V8- On the other hand, both the calculations and the simulations 
show a large peak in g(uj) around uj — 0.75. This peak is not present for 
this density (9%) of randomly placed weak springs, as found by^ and also 
in our own calculations and simulations for uncorrelated weak springs. It is 
very reminiscent of the peak in the density of low frequency vibrations found 
in glasses, where it is known as the Boson peakP The difference in the heights 
of the maxima between the calculations and simulations are associated with 
the smoothing that has to be performed for the latter, and the areas under 
the peaks, which correspond to the total extra density of states in this region, 
are very similar. The slight difference in the positions of the peak may also 
be associated with this smoothing. However, it might just be a genuine effect 
arising from interactions between the defect clusters which are ignored in the 
AVTA but could occur in the simulated systems (where noise places a lower 
limit to the defect concentrations that can be studied). Incidentally, in the 
calculated density of states, there are slight problems at the frequencies of the 
Van- Hove singularities, as can be seen in figure 1, because of our use of equation 
(17), since the derivative of Gq(z) is discontinuous at these points. 

The main conclusions from this paper are that the presence of defect clusters 
involving correlated disorder extending over several sites can have an apprecia- 
ble effect on the density of states. The average trace approximation (AVTA) 
presented in this paper provides a simple (at least in principle) and accurate 
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method of calculating the density of states in such systems for low concentra- 
tions of such clusters. 
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Caption for figure. 

Fig.l (Color on line) The density of states g(u>) as a function of the frequency 
u> for a system containing a concentration 1.5% of defects with six weak springs 
all having K\/Kq = 0.1. The continuous line shows g(w) for the system without 
defects, the dashed line the density calculated using equation (17) and the dotted 
line that obtained from the simulations. 
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